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Abstract 

In this paper we calculate at high-precision the Laurent expansions in e = (4 — D)/2 of the 
17 master integrals which appeared in the analytical calculation of 3-loop QED contribution to 
the electron g-2, using difference and differential equations. The coefficients of the expansions 
so obtained are in perfect agreement with all the analytical expressions already known. The 
values of coefficients not previously known will be used in the high-precision calculation of the 
4- loop QED contribution to the electron g-2. 
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The g-2 of the electron is probably one of the most precise test of QED. Over the years, the con- 
tinuous improvements in the precision of experimental determinations have demanded correspondent 
improvements of the theoretical predictions. Let us recall the current status of calculations of QED 
contribution to the electron g-2. One- and two- loop contributions are known in closed analytical form 
for a long time. The calculation of three-loop contribution in closed analytical form was completed 
more recently [1], after many years of hard work. Four-loop contribution is known at present only 
numerically [2], with a precision of about 2.5%, obtained using Monte-Carlo integration methods. At 
present this precision is adequate for comparison with current experimental determinations. Anyway, 
there is the need of a new independent high-precision calculation, in order to cross-check the current 
numerical value and to improve considerably its precision, in view of future experiments. 

A technique which turned out to be useful in g-2 calculations is integration- by-parts [3, 4]. The 
contribution to g-2 of a graph, expressed in the form of a combination of many different integrals 
with different powers in the numerator and in the denominator, is reduced to a linear combination of 
a small set of "master integrals" by using identities obtained by integrating-by-parts in D-dimension 
space-time. The master integrals must be calculated analytically or numerically in the limit D — > 4 
by some method. 

In Ref. [1] this technique was applied to the analytical calculation of the contribution to the elec- 
tron g-2 of the last family of three-loop graphs still not known analytically, the so-called triple-cross 
graphs; the contributions of all the other three-loop graphs were already obtained in analytical form 
by other methods. The contribution of triple-cross graphs to the g-2 was reduced to a combination 
of 18 master integrals, called I\, J 2 , . . . , I\&. The Laurent expansions in e = (4 — D)/2 of the master 
integrals were calculated in analytical form, by direct calculation (the simpler ones) or by using 
identities which relate the coefficients of expansions to values of integrals in 4 dimensions already 
known from previous work[5, 6, 7, 8]. Subsequently, in Ref. [9], we found that the QED contribution 
of all three-loop graphs can be reduced to a linear combination of the same master integrals, which 
are therefore the only master integrals needed in the three-loop calculation. We also found that 
one master integral, In, is a linear combination of the integrals 1 14 and Jig so that only 17 master 
integrals are needed. The topology of the 17 master integrals is shown in Fig.l. 

We plan to use the integration-by-parts method in the calculation of the 4-loop contribution. We 
expect about 300 ~ 400 master integrals; an analytical calculation seems to be out of reach, so that 
alternatively we consider a high-precision numerical calculation. Of course, we begin by calculating 
the simplest master integrals, those which factorize in products of master integrals with fewer loops. 
At 3-loop level there are two master integrals, /15 and lie (see Fig.l), which factorize in one 2-loop 
master integral multiplied by one one-loop tadpole. At 4-loop, analogously, there are 17 4-loop master 
integrals which factorize in the product of one 3-loop master integral Ij and the one-loop tadpole. 
The expansion in e of the one-loop tadpole contains 1/e, so that we need one additional term in the 
expansions of the 3-loop master integrals in order to obtain 4-loop master integrals expanded at the 
same level required in the 3-loop calculation. Unfortunately the number of terms of the expansions 
of Ij of Ref. [1] just suffices to obtain the 3-loop result. 

Therefore, for this reason, in this paper we have calculated at high precision (30 digits) deeper 
expansions in e of the 17 master integrals If 

hC~ l = -.74727427517503872293442889217- 5.3830601530185743669980288649e 

- 26.213893908510814642139484925e 2 - 105.674804179662879301749904804e 3 

- 385.27172279141197066520099021e 4 - 1318.9119175947583164248019535e 5 

- 4337.082974360622762021390292e 6 + 0(e 7 ) , (1) 
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l 2 C- x = 2.4041 13806319188570799476323026- 1 - 5.33950114034530650865549677321 
+ 25.5217203762232230112382243421e - 73.5271290561422487313093251595e 2 
+ 249.044154178834576505276496985e 3 - 753.473335692793671251435382464e 4 
+ 2346.37248004908075896840609088e 5 - 7109. 19379761589640023765834979c 6 + 0(e 7 ) , (2) 

I Z C' 1 = 0.333333333333333333333333333333e~ 3 + 2.33333333333333333333333333333e~ 2 
+ 10.33333333333333333333333333336" 1 + 19.7427119912535493212749714262 
+ 88.8774427983382280511485935554e + 95.6785746389971388003677314582e 2 
+ 626.528599797078129746832638577e 3 + 210.711024315707077431312202546e 4 
+ 4305.95682691978475295985357628e 5 - 1994.41925210304147109855938156e 6 + 0(e 7 ) , (3) 

hC' 1 = 2.404113806319188570799476323026" 1 - 1.88228069596232589686771921947 
+ 24.9738222767976352300663548079e - 40.24308322651 14989313433432375e 2 
+ 216.453265450196273372723759103e 3 - 490.48300355255485916860004085e 4 
+ 1877. 77603630120111667601080315c 5 - 4980.53077825841233320374729215e 6 + 0(e 7 ) , (4) 

hC- 1 = 0.166666666666666666666666666667e~ 3 + 1.5e~ 2 

+ 5.876798532970213793721836333376" 1 + 10.2027383130242843587260466742 

+ 49.4062940751534383431822667245e + 38.7587897454906373636992141483e 2 

+ 355.54075314113673706462779833e 3 - 12.8313984708255340018628252229e 4 

+ 2544.39228074929353310098747381e 5 - 2170.06912961044909561870930877e 6 + 0(e 7 ) , (5) 

/gC- 1 = 0.333333333333333333333333333333e~ 3 + 2.33333333333333333333333333333e~ 2 
+ 10.3333333333333333333333333333e _1 + 29.7659868661137435310832928686 
+ 102.438098206123109827196396063e + 229.885864095858993939826105646e 2 
+ 762.395106594671638134114238656e 3 + 1423.01797639907075045462241083e 4 
+ 5052.27961041516425974606724499e 5 + 7498.83181107570646929366512712e 6 + 0(e 7 ) , (6) 

I 7 C- 1 = 0.166666666666666666666666666667e~ 3 + 1.5e~ 2 

+ 5.876798532970213793721836333376- 1 + 18.0613667204622021855023607236 

+ 55.4913080622726427983915540249e + 134.432560982177142077741730996e 2 

+ 394.304047315939025666187053151e 3 + 824.935548393665344225617774097e 4 

+ 2505.10730543388776893266213278e 5 + 4469.58551067612956096036609952e 6 + 0(e 7 ) , (7) 

hC~ l = -e~ 3 - 5.33333333333333333333333333333e~ 2 - 16e _1 

- 43.9148312632524344127591663433 - 154.918028662961860140720725441e 

- 374.094185333355168755155126113e 2 - 1436.67271253544219472250003758e 3 

- 3281.94043631920573984798601407e 4 - 13100.0652782910379543889228475e 5 

- 29197. 1252330882109400594603379e 6 + O(e 7 ) , (8) 
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IgC- 1 = -0.666666666666666666666666666667e~ 3 - 3.33333333333333333333333333333e- 2 

- 11.956534800363119539611497e- 1 - 44.5995196208454844578584038609 

- 125.546733075809915579874065329e - 442.029796112593386030203112952e 2 

- 1165.46537028123770794611984172e 3 - 4078.5492471997365803005187143e 4 

- 10533.3774099326183047654759014e 5 - 36961.4052223824089107170672131c 6 + 0(e 7 ) , (9) 



I^C- 1 = -0.333333333333333333333333333333e~ 3 - 1.66666666666666666666666666667e~ 2 

- 10.5797362673929057458896606666e~ 1 - 30.1135700965916539174115430661 

- 136.273953266205829743326933622e - 333.735442456113204996932778211e 2 

- 1376.89949509222838398713383632e 3 - 3194.567378327629536037718821e 4 

- 12862.3364927028759414511191045e 5 - 29234.6901535778480407834346202e 6 + 0(e 7 ) , 

(10) 



l X2 C- x = e~ 3 + 3.5e~ 2 + 7.027777777777777777777777777786- 1 

+ 11.5787037037037037037037037037 - 24.3219708028082785734403310608e 

- 121.961009709497309747947635538e 2 - 819.989946014746295714157758046e 3 

- 2404.86316062636537480556898164e 4 - 10421.8794385530161287164813459c 5 

- 27852.0785333053790606284709671e 6 + 0(e 7 ) , (11) 



I lz C~ l = 2e~ 3 + 7.66666666666666666666666666667e- 2 + 17.56- 1 

+ 22.9166666666666666666666666667 + 21.2517910512915199882568913631e 

- 184.230005105298483421085631284e 2 - 661.110586153353833254708841052e 3 

- 3685.05477938169880569422657784e 4 - 10050.975403938380587540545542c 5 

- 42319.9745464132618093606637193e 6 + 0(e 7 ) , (12) 



I^C' 1 = 0.333333333333333333333333333333e~ 3 + 1.16666666666666666666666666667e~ 2 
+ 2.08333333333333333333333333333e~ 1 + 2.99715174175891809439930176403 

- 21. 7478477083774449677944881844c - 84.9528451829928094897709758796e 2 

- 520.302648810090620535897692307e 3 - 1504.03113911470006499851 17616e 4 

- 6482.74483328181932000540551865e 5 - 17315.8330540574013027834356625e 6 + 0(e 7 ) , 

(13) 



I^C- 1 = 1.5e~ 3 + 5.75e~ 2 + 13.1256- 1 + 30.3469725347858114917793213332 

+ 59.9061795885052938207577452894e + 128.784074758885330304917817688e 2 

+ 247.030559404244292769216590768e 3 + 522.536074129783764428998974715e 4 

+ 995.526049616742468296298492052e 5 + 2097.54543056218400478046531901e 6 + 0(e 7 ) , 

(14) 
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I W C^ = 0.5e~ 3 + 1.75e~ 2 + 6.414868133696452872944830333296- 1 

+ 16.0102660805759621969058588126 + 42.0560659585935725800136421336e 
+ 91.5580001461564241262070311808e 2 + 211.940358797491166303979686218e 3 
+ 432.923773655861990482605036606e 4 + 946.511285468009699310821028163e 5 

+ 1874.61717426917014415736822215e 6 + 0(e 7 ) , (15) 

InC' 1 = -0. 166666666666666666666666666667c" 2 - 0.972222222222222222222222222222c" 1 

- 5.87783109665941583590779329626 - 23.4899848806977618982607310895e 

- 97.5806595141990146514480659561e 2 - 330.869679035399716046952939963e 3 

- 1181.19322772402976499602304607e 4 - 3694.74722374791467371175473492e 5 

- 12281.2292376715157620002229398c 6 + 0(e 7 ) , (16) 



I^C- 1 = -e" 3 - 3e~ 2 - 6e _1 - 10 - 15e - 21e 2 - 28e 3 - 36e 4 - 45e 5 - 55e 6 + 0(e 7 ) . 

(17) 

The integrals Ij are defined as[l] 

I j=(^) I d D hd D k 2 d D k 3 if , P l=p .k 2 , P j = liij>2, 



3 \n D ~ 2 J J " "* Qj ' 

Qj contains the product of the denominators of the corresponding j-th graph of Fig.l. The normal- 
ization factor is C = (7rT(l + e)) 3 . The number of terms of the expansions of Eqs.(l)-(17) suffices 
for the use in g-2 calculations at four and even more loops (actually, we have calculated much deeper 
expansions, not completely listed here for lack of space). 

Eqs.(l)-(17) are to be compared with the corresponding analytical results of Ref.[l]. Due to an 
unfortunate misprint (see [9]), in Ref. [1] the terms containing the constants C±, C 2 are missing in the 
r.h.s. of the integrals I 2 , h, h, I5, I7 and In; all the results which follow are however correct, as 
Ci, C 2 cancel out systematically in the final results of [1]. For the ease of the reader we list here the 
correct expressions of the integrals I 2 , J 3 , . . . , I 7 (already appeared in [9]); the remaining integrals 
are listed in Ref.[l]: 

, 1 n C(3) 13 4 1 2 , , f 385 ^ 85 2, 



I2C- 1 = - ^ - 3^ + 10C(3) + el — C(5) - -tt 2 C(3) 



- — 7T 4 - 82C(3) - 47r 2 ln 2 + 16tt 2 - 2Ci + QC 2 J + 0(e 2 ) , (18) 
15 J 



4((3)-8A2+^ + ? + 4C 2 ] fO(r). (19) 
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rr -i 9 C(3) 7 4 1 2 /385 85 , 7 4 

AC? = 2— --vr +2C(3) + - 7 r + e^— C (5) - -n C (3) - -tt 



82C(3) - 4tt 2 1ii 2 + 16tt 2 - 2d + 4C 2 ) + 0(e 2 ) , (20) 



rr i 1 3 1 1 55 \ 4 4 14 7 , 95 



+ c( -|tt 4 - 44C(3) - |vr 2 + ^ + 2d ) + 0( e 2 ) , (21) 



i 1 7 31 4 4 2.,„, 1 2 103 /45., . 7 2 . /0 , 11 4 

J 6 <? = — ? + — ? + VT + -((3) + -7T 2 + + e — C(5) - -7T 2 C(3) + TT 

3e 3 3e 2 3e 45 3 SV ; 3 3 \ 2 2 45 

+ y C(3) - 4vr 2 ln 2 + ^r 2 + ^ + 2d) + 0(e 2 ) , (22) 

T „ i 1 3 1/ 1 2 55\ 1 4 8 >/ x o 95 /45 w , 17 2>/ , 

- V - 50((3) - 47T 2 ln 2 + i» 2 + ^ + 2C 2 ) + 0(e 2 ) , (23) 
9 o o y 

where = X^iV nP > °4 = Z)nLi an< ^ ^ 2 are ^ ne coefficients of the e term in the 
expansion of the integrals 

-% \ 3 f d D k 1 d D k 2 d D k 3 (p.jfci) / 2 4 _ _ , 2 , 



d D hd D k 2 d D k 3 ( P .k 3 ) ri ( l 4 1 , , ., 2 



m > - i^j j - c 1-30- - r + ^ + + °n 



£>i = (p - A^) 2 + l - i0, D 2 = (p - fci - A; 2 ) 2 + 1 - i0, L> 3 = (p - &i - A; 2 - A; 3 ) 2 + 1 - iO, 
D 4 = (p-k 2 - k 3 ) 2 + 1 - i0, D 5 = (p - k 3 ) 2 + 1 - i0, D 6 = k\ - i0, D 7 = k 2 2 - i0, D 8 = k\ - i0, 
p 2 = -1. 

The analytical values of C 2 and d were calculated, respectively, in Ref. [9] and Ref.[10] 
d = -^C(5) + gvr 2 C(3) - ^4 + 18C(3) + 2?r 2 ln2 _ 37r 2 

Ci = -f C(5) + |vr 2 C(3) - ^ + C (3) + 2vr 2 In 2 - £ . 

Eqs.(l)-(17) agree perfectly with the analytical expressions of Eqs.(18)-(23) and of Refs.[l, 9], 
and with the deeper expansions of ii [10], J 15 [ll], J 16 and In- 

Now we sketch the method used for obtaining Eqs.(l)-(17). We have used in this calculation the 
new method of calculation of master integrals based on solution of difference equations in exponents 
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developed and described in detail in Refs.[12, 13]. This method consists in the construction and nu- 
merical solution of systems of difference equations between the master integrals Ij (with polynomials 
in D and n as coefficients), seen as functions of the exponent n of one denominator. The difference 
equation for a given master integral contains in the r.h.s. only master integrals with less denomina- 
tors, which are simpler; therefore the equations of the system can be resolved one at once, beginning 
with that corresponding to the simplest master integral and ending with that corresponding to the 
most complex. Suitably boundary conditions at n — > oo must be provided for a proper solution of 
the difference equations. This implies the calculation of integrals with one loop less, which can be 
obtained by solving other systems of difference equations. The amount of calculations needed to work 
out and solve the systems of difference equations is rather high, so that the practical application of 
the method relies on the use of an automatic tool, the program SYS described in Ref. [12]. 

Because of some actual limitations of the program SYS, we have found it convenient to give a mass 
A 7^ to the photon lines, so that the master integrals become functions Ij(X). Then, we have used 
the above described method to calculate the value of master integrals (actually, these values 

were already calculated using this method in Refs.[14]; here we have repeated the same calculation 
increasing the numerical precision and the order of expansions in e). 

Subsequently, always by means of the program SYS [15], we have built a system of differential 
equations in A for Ij(X) (the approach based on the use of differential equations in masses was 
introduced in Ref. [16]), and we have integrated it, using as initial conditions the values (and, when 
needed, the derivatives w.r.t. A) at A = 1, obtaining the values 7,-(0), that is, Eqs.(l)-(17). 

We note that the systems of difference and differential equations for A ^ are much more 
complicated (more master integrals, higher order equations, higher degree coefficients) than the 
system of difference equations for the original master integrals with zero masses lj(0), so that the 
solution of the latter would be preferable for calculating 4-loop master integrals which do not factorize. 
We plan to overcome this limitation in near future. 
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Figure Captions 

Figure 1: Topologies of the master integrals. Dotted lines are massless. 



(S. Laporta, High-precision e-expansions of three-loop master integrals contributing to the elec- 
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